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1. This paper treats of the transverse vibrations of a circular disk of 
uniform thickness rotating about its axis with constant angular velocity. 
The problem has a practical bearing, as throwing light on the occasional 
failure of turbine disks. In an impulse turbine blades are fitted to the rim 
of a thin disk, and steam is admitted to them from nozzles which are usually 
arranged symmetrically, but not always continuously, round the periphery. 
The failures appear to be sometimes due to the blades coming in contact with 
the adjacent parts. This points to transverse vibrations, and there can be 
little doubt that the phenomenon is one of resonance between the periodic 
forces exerted by the steam jets and the periods of free vibration of the disk.f 

For a thorough investigation it would be necessary to ascertain, to a 
reasonable degree of approximation, the free periods of asymmetrical disk 
of any given profile. This may form the theme of a subsequent investiga- 
tion by one of the authors. In the meantime the study of the simpler 
case now discussed illustrates the physical aspects of the problem, and has 
suggested an important simplification in method, 

2. The special difficulty of the problem arises from the fact that centrifugal 
force and flexural rigidity are both operative, and that the normal modes 
of vibration depend on their relative importance. The vibrations of a non- 
rotating disk were fully discussed by KirehhofT in two well-known papers.J 
In the other extreme, where the disk is so thin and the rotation so rapid that 
the flexural forces are negligible, the problem admits of a very simple 
solution, as follows : — 

When an elastic disk rotates in its own plane, in relative equilibrium, the 
principal tensions are along and perpendicular to the radius vector. If a be 

* My attention was -drawn to this question by Prof. Gerald Stoney, who showed me 
some pretty experiments on a rotating india-rubber membrane. This led to various 
calculations, I found subsequently that Mr. Southwell had become interested in the 
same problem, and a comparison showed that the methods and the results were closely 
parallel. This paper is accordingly a joint production, but it should be stated that the 
important remark that the formula p 2 =p x 2 +p 2 2 °f § 4 g ives not merely an approximation 
but a lower limit to the frequency is due to Mr. Southwell.-— H. L. 

t A paper by K. Baumann, read before the Institution of Electrical Engineers on 
March 15, 1921, gives examples of failures which have been traced to this source. 

J "Ueber die Schwingungen einer elastischen Scheibe," ' CrehVs Journal, 5 vol. 40 
(1850), and « Pogg. Ann., 5 vol. 81 (1850). 
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the radius, p the density, and co the angular velocity of rotation, their values 
at a distance r from the centre are* 

P = A(a 2 -r 2 )pw 2 , Q = (A^-Br 8 )/)©*, (1) 

respectively, where 

A = i(3 + <r), B = i(l + 3<r) t (2) 

<t denoting Poisson's ratio. These values are obviously unaffected by a small 
displacement (to) normal to the plane of the disk. 

The equation of transverse motion of an element of area whose polar 
co-ordinates are (r, 6) is, under the condition stated, 

' rt *ff-s( w »*)* + »( , *a) M - <3) 

01 P dt 2 " r 9r \ 3r / + r 2 30 2 ' W 

Assuming that w varies as cos s0 cos pit, where 5 is integral, we find 

A^-^f^+l |f)-2Ar ^ + (P^ + ^B-^f) w = 0. (5) 

\ar J r cr / or \(o u r^ j 

This is readily integrated by a series. Assuming 



rp \fC 



w = tC k [L), (6) 



we have 



C k+2 _ /g (A + 2) A -s 2 B -pS/v 2 
V k {(k+2f-s 2 }A ' 



(7) 



or, if we write 



£L + s 2 £ _ m ( m + 2)A, (8) 



2 

C^+2 _ (Jc—in)(Jc-hm + 2) 



C k "(fe + 2-*)(* + 2 + «y (9) 

The solution which is finite for r = begins with (r/a) s . Hence, inserting 
the factors temporarily omitted, 

W 1 2(2s + 2) a s 

2.4 (2s + 2) (2s + 4) a 4 J n \ / 

The series in { } is hypergeometric. In Gauss's notation, 

A) F \ ' A ' ' a> 

where a = J(s— m), ft = -| (s -f- m -f 2), 7 = 5 + 1. (12) 

* Love, * Elasticity,' § 102 (6). 

x 2 



w?== ^(;;) **> & -v> Ti) COSS^COSpi^, (11) 
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The expansion is obviously convergent for r < a, but since 

7 = « + A (13) 

it. becomes, by a known theorem, logarithmically infinite at the edge r = a, 
unless it terminates.* Hence we must have 

m = s + 2?^ ) (14) 

where ?i is a positive integer. Thus 

\al 1 1(8+1) a 2 

4-i L^ ^^L± ! / — — ... > COSS0COSPi£ (15) 

1.2(s+l)(s + 2) a 4 J r v y 

Arbitrary constants may of course be added to 8 and £. 

The mode of vibration thus found is characterised by n nodal circles and 
s equidistant nodal diameters. The corresponding frequency is deter- 
mined by 

^ = (s + 2,n)(s + 2ri+2) A-^B, (16) 



where A and B have the values (2). That this ratio should be independent 
of the diameter and thickness of the disk might have been anticipated from 
the theory of dimensions. 

The modes for which s = involve displacement of the centre of the disk,f 
and are therefore not relevant to the case of a disk mounted on a shaft. 
The latter remark applies also to the mode where s = 1, n = ; the disk then 
remains plane in form, but is slightly untrue to the axis. The most inter- 
esting modes, and the ones most easily realised experimentally, are those in 
which there are no nodal circles. Putting n = 0, we have 

w = C( — ) coss0cosjpi£, (17) 



Pi 



a i 

2 



with ¥\ = £(l-<r)tf» + J(3 + <r)s. (18) 



When s = 2, we have two nodal diameters at right angles, and 

^ = i(5-^ (19) 

Por an india-rubber membrane we may put <r = |, practically, whence 
pfw = 1*5. 

* Of. Bayleigh, 'Theory of Sound/ §336. 

t It may be noticed that if we put rja— sin|<£ the formula (15) may in this case be 
written 

W =. CP W (COS </)) COS jDj£, 

where ~P n is the zonal harmonic of order n. 
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Writing (16) in the form 



J?i 2 = Xo) 2 , 
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(20) 



the values of X for various modes are given in the following table, on the 
hypothesis a = 0*3, which corresponds fairly to the case of a steel disk. The 
mode for which n = 0, s = 2, is the most important from the practical point 
of view, since it is the gravest type of free vibration if we exclude (for the 
reason already given) the case of n = 0, s = 1. 



n. 


* = 1. 


*»2. 


5 = 3. 





1 


2 35 


4-05 


1 


5 -95 


8-95 


12*3 


2 


14 2 


18-85 


23*85 


3 


25-75 


32'05 


38-7 



It is to be observed, throughout this paper, that the vibrations considered 
are vibrations relative to the undisturbed disk ; in particular, the nodal lines 
in the modes investigated are in reality lines of rest in the disk, and are 
accordingly carried round by the rotation. The results may, however, be 
stated in a different form. By superposition of two normal modes of the 
same type, with a suitable difference of phase, we can construct wave-forms 
which travel round the disk with the angular velocity ±pifs. The preceding 
results would, in fact, have been unaltered if we had assumed that w varies 
as cos(pi£ + s#) instead of cos£>i£coss#. The corresponding angular wave- 
velocities in space are oy+pi/s. Thus, in the case n = 0, s = 1, where p\ = o>, 
by (16), the wave-velocities are 2a> and 0. In the latter alternative, the 
membrane merely rotates in a fixed plane slightly inclined to that which 
we have taken as the plane of reference. 

3. The preceding may be regarded as results which are approximated to as 
the speed of rotation is increased. In the other extreme, where the influence 
of the rotation is negligible, and the flexural rigidity predominant, we may 
have recourse to the investigation of Kirchhoff. The results are of the form 

EA 2 



iv = a* 



pa* 



(21) 



where E is Young's modulus, h is the half-thickness, and //, is a numerical 
coefficient depending on the particular mode of vibration and (to some 
extent) on the value of o\ The results are actually presented by Kirchhoff 
somewhat differently. He gives tables of frequencies of various modes 
relative to that of the gravest mode of all as a standard, first on Poisson's 
hypothesis of <r = -J, and secondly for a = £. The gravest mode in question 
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has two nodal diameters and no nodal circle. The frequency for this is 
given by Kirchhoff in semi- vibrations per second, viz., 



and 



p2 
7T 

7T 



V(^> 1-04604, 
V(^)xl-02357, 



(22) 



(23) 



for a = \ and a = ^, respectively. 

An independent calculation for the case of <r = 0*3 gives the following 
values of the coefficient jul in (21). The number of nodal diameters is 
denoted, as before, by s, and that of nodal circles by n* 



n. 


s = 0. 


s = 1. 


s = 2. 


s = 3. 









10-61 


56-87 


1 


29*67 


154-3 


443-9 


1040 


2 


577 


1222 


2624 


4285 


3 


2983 


5266 


8574 


13287 



4. When the centrifugal tensions and the flexural stiffness are 
effective, the difficulties of an exact solution are much increased, 
equation of motion is 



both 
The 



c/w 



p dfi ~ • 



■r dr 



Vr 



JW 



dr 



+ 



Q &w 



Wi> 



dd 2 3(1 -a 2 ) 



V%; 



(24) 



and the boundary conditions are as given in Kayleigh's " Theory of Sound," 
§ 215. If we assume that w varies as cos s6 cos pt, as before, the equation 
can be integrated by a series, but the relation between successive coefficients 
consists of three terms, and the deduction of general results would not appear 
to be easy. 

Approximations which are adequate from the practical point of view can, 
however, be obtained by Kayleigh's method of an assumed type,f at all 
events, for the more important modes. The procedure is equivalent to this, 
that we calculate the potential and kinetic energies on the basis of an 
assumed form oiw, and determine the frequency so that the mean energies 
of the two kinds shall be equal. Since the frequency thus determined is 
stationary for slight variations of the type, a good approximation can be 
found in this way if the type be suitably chosen. It is known, moreover, 
that in the case of the gravest mode of all, the frequency obtained will be an 

* The authors are indebted to Miss B. S. Gough for assistance in the rather laborious 
calculations. 

t < Theory of Sound/ §§88, 89. 
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upper limit to the true value. In the present problem, this statement 
obviously applies also to the gravest mode having any assigned number s of 
nodal diameters. 

The potential energy of the centrifugal forces is 

Vl = h \\i T &f +Q $&f} rd0dr> (25) 

where P and Q have the values (1), and the integration extends over the area 
of the disk. 

The potential energy of flexure is, in Cartesian co-ordinates,* 

Y 2 = - q/ti 2\ {(yhv) 2 + 2{l~<r)(io^ 2 --w xx Wyy)} dxdy. (26) 

To transform to polar co-ordinates we have 

V 2 W _ ^ + l^f + I§!!f, (27) 

dr 2 r dr r 2 3# 2 ' 



1 d 2 w 1 dw\ 2 d 2 w 1 1 a 2 w , 1 div\ 

V 



2 — / x u w x vw \ u w / x u w i 

ay yy — ^ g^g^ ^ ^g J ^ 2 \ 2 g^ 2 ^ g ? . 



The kinetic energy is 

T = P h[[(^)*rd0dr. (29) 

If, then, we assume 

vj = <f> (r, 6) cospt, (30) 

the method referred to leads to a result which may be written 

V 2 _ Vi W + v 2 W (m 

p __ — , {i i } 

where T (<j>) differs from T in that dw/dt is replaced by <£. 

In estimating in this way the frequency of the gravest mode for a given 
value of 5, it is sometimes convenient to adopt for <f> a form containing a 
variable parameter, and to adjust this afterwards so as to make the fraction 
in (31) a minimum. This reduces the upper limit referred to. This plan is 
adopted, and illustrated by numerical examples, in § 5 below. The work is 
somewhat intricate, but the formula (31) suggests a simpler process, which is 
probably adequate for practical purposes, having regard to the fact that the 
elastic constants E and <r for a given disk are not likely to be known to a 
high degree of accuracy. 

It is evident that if precisely the right form be given to <j> the formula (31) 
will be exact. Now, on account of the stationary property referred to, the 
true value of <f> would give an approximate value for pi 2 , viz., 

*• = w (32) 

* ' Theory of Sound,' §214. 
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In the same way the true value of </> would give an approximate value for pi, 
viz., 

* - w>- <33> 

Moreover, by Kayleigh's theorem, these estimates will in the case of the 
gravest mode having an assigned number of nodal diameters be too high. 
The second member of the equation (31), when <£ has the correct value, is 
therefore greater than p x 2 +p 2 2 . Hence in the case supposed the formula 

p 2 = pi 2 +P2 2 , (34) 

where p x and p 2 are the values of p on the two extreme suppositions of infinite 
thinness, and infinitely slow rotation, respectively, gives an approximation 
which is of the nature of a lower limit to the true value.* 

The property here inferred may be generalised as follows. If the restoring 
forces which control the vibrations of an elastic system can be separated into 
two or more groups which affect the potential energy independently, and if 
the gravest frequency of vibration be found on the assumption that each 
group acts separately (the inertia being unchanged), then the sum of the 
squares of the frequencies thus found is less than the square of the greatest 
frequency which can occur when all the groups act simultaneously ; i.e., it is 
less than the square of the gravest frequency natural to the system. 

As a numerical example, suppose that 

E = 2 x 10 12 , p = 7-8, a = 60, h = 1, © = IOOtt, 

in O.G.S. units. In the gravest mode with two nodal diameters we have, on 
the hypothesis that cr = \, 

%£ == 2-375 x 10 4 , fi? = 2-16485 x 10 4 , 

from (19) and (22). The formula (34) gives 

P = ViPl±£Jl = 213-07, (35) 

IT IT 

as a lower limit to the frequency. 

5. We proceed to illustrate the application of Kayleigh's method to the 
problem. 

We assume 

io = 1 4-/3 -s ) ( - ) cos s8 aoaptf (36) 

*'' See footnote on p. 272. 
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where /3 is a constant to be afterwards chosen so as to make the inferred 
frequency a minimum. We find 

, r , 2 27 rs(s + 2)A— s 2 B , 0O s(s + 2)A-s 2 B 

Vl = * •*P t * a * h I s+1 + 2 # ^2 

+ /S (« + 2)(.+ 3) J 008 ^' (37) 

V 2 = 7r ^ 2 3 )ffi2 [8(s+l) j 3 2 +2(l-a){^(a-l) + 2/3 S (^-l) 

+ /3 2 (s + l)(s 2 -2)}.]cosV, (38) 

T = i 7 rM« 2 (-^ T + ^ + -^)^ 2 sin 2 ^. (39) 

\s+JL s + ^ s + o/ 

Equating the mean value of Vi + V 2 to that of T, or expressing that the 

total energy is constant, we have 

2 = L + 2M/3 + N/3 2 (m 

1 l + 2m/3 + n/&' 

where L, M, 1ST, I, m, n have the values given below. "We have to choose /3 so 
as to make this fraction a minimum. 

The stationary values of p 2 are given by the quadratic 

(Z»-w 2 )|?*-(nff+.»L-2wM)y + LN--M 3 .= 0, (41) 

and the corresponding values of y3 by either of the equations 

S - M 7" m -P 2 = ^~¥ m (42) 

wp 2 — N mp 2 — M 
We write for shortness 

C» = ^L, (43) 

per 

so that f (like g>) is the reciprocal of a time. The coefficients in (40) have 

then the values 

L = S (s + 2)A-s 2 B w2 + 4s 2 (s-l) ?2; (44) 

s + 1. 3(1 + 0") 

M = g ( g + 2)A- S 2 B &)2 + 4ay-l) f2; (45) 

s -+- 2 3 (1 + a") 

„ _ , ( s 3 + 4a 2 +4g+4)A-g 2 (s + 2)B . 2 16(s+l)+4(l-<r)(g + l)(s 2 -2) „ 

(s + 2)(s + 3) W+ 3(l-a 2 ) fe ' 

(46) 

/ = — m = ^, == (47) 

s+1' s + 2' s + 3 v y 

We have specially in view the case of s — 2. Adopting the value J for cr, 
we have 

A = I B = l? (48) 
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Hence 

L -2i m+ To^ M -32 M+ ^ ? ' ^-160 W+ U ? ' (4y) 

l = \, m =\, » = 1. (50) 

As a partial verification of the method we note that when to = the 
quadratic (41) reduces to 

15( / ^Y~19712|^ + 212992 = 0, (51) 



the lower root of which is 



This makes 



f 



*L = 10*8955. 
f 2 



£ = 1*05069 £ (52) 

IT 

This is higher than the correct value (22), as was inevitable, but the error 

is less than J per cent. 

To illustrate the case where the effects of centrifugal force and flexural 

rigidity are comparable, take the numerical data given at the end of § 4. 

These make 

co 2 = 9*86960 x 10 4 , £ 2 = 1*97847 x 10*, (53) 

L = 1*62549 x 10 5 , M = 1*85223 x 10 5 , N = 519182 x 10 5 . (54) 

The quadratic (41) becomes 

0-00417p 4 -ri2959 xl0 5 p 2 + 5*00850 x 10 10 = 0. (55) 

The lower root is 

f = 4*50890 x 10 5 , (56) 

whence, for comparison with (35), 

£ = 213-74. (57) 



7T 



This is probably very close to the true value. In any case the upper and 
lower limits which we have obtained differ only by about \ per cent. For 
practical purposes the approximation (35) would probably be sufficient in all 
cases, for reasons given above. It is, of course, obtained with much less 
labour than (57). 

The value of fi corresponding to (56) is —0*169. 



